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Abstract 

We study interaction of low mass dark matter within beam dump experiments. In particular 
we study the dipolar dark matter model which assumes that the dark matter couples to Standard 
Model particles via its electric or magnetic dipole moment. We analyse the constraints on this 
model in the context of a particular beam dump experiment E613 conducted in the Fermilab. We 
find that dark matter mass in the range of 1 — 10 GeV with a magnetic dipole moment between 
(0.33 — 1.5) X pb and a electric dipole moment between (0.5 — 3) x 10“^^ e-cm. We compare 
the bounds from other experimental data, such as helioseismological data and direct detection 
experiments. 
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I. INTRODUCTION 


Evidence from galactic rotation curves as well as bullet cluster collisions have strongly 
suggested the existence of dark matter (DM). However the nature of DM is unknown. In 
theory there exist many candidates for DM, the most popular of them being the case of 
weakly interacting massive particle (WIMP). The WIMP scenario is realized in many models, 
with neutralino in supersymmetric theories being the most well studied of them. However, 
there also exists the possibility of DM interacting electromagnetically with Standard Model 
(SM) particles [Tl-fTn]. We refer to this model of DM interactions as dipolar DM model. 

The dipolar DM model assumes that dark matter couples to photons through loops to 
give rise to electric and magnetic dipole moments. Here we study the case in which DM 
particles are Dirac fermions. The effective Lagrangian for such a DM particle interacting 
with an electromagnetic held (-T/ii^) through its electric dipole moment (D) and magnetic 
dipole moment (/r) is 


= + ( 1 ) 

Here we consider this particular scenario in the context of DM detection in so called beam 
dump experiments also known as hxed target experiments. Complimenting the searches 
from direct and indirect detection experiments, beam dump experiments feature a high 
energy beam incident on a hxed target. The DM particles produced from this collision 
are then detected in a suitable detector. In this work we particularly focus on the E613 
experiment at Fermilab in which a 400 GeV proton beam is incident on a tungsten target 
with the resultant DM produced from the annihilation of beam proton and target proton 
being detected in a lead detector after passing through iron shielding. The advantage over 
direct detection in such experiments is that of higher luminosity. But the reach of such 
experiments and in particular of E613 is restricted to low mass DM due to the kinematics of 
hxed target experiments. Future experiments like the new hxed target facility proposed at 
the CERN SPS called SHIP (Search for Hidden Particles) [Hj can explore this possibility. 
The possibility of DM detection at beam dump experiments has been studied in the case 
of light mediators in the dark sector mediating DM interactions with SM particles [HI OS] 
and also in the case of Z' as the mediating particle HI. More recently a similar model with 
dark vector boson was studied in the context of E613 Ha. We follow the approach in Ha 
and study the constraints on the dipolar DM model from E613. 

The paper is organized as follows : In section 2 we describe the method for calculating the 
production cross section of dipolar DM when the 400 GeV proton beam strikes the tungsten 
target. In section 3 we describe the deep inelastic scattering that takes place between the DM 
and the lead target nuclei. In section 4 we give the results from analysing the dipolar DM 
in the case of E613 and compare constraints from other experiments. Finally we conclude 
in section 5. 


II. DIPOLAR DM PRODUCTION IN E613 

Here we follow the procedure outlined in [13]. In a beam dump experiment like E613, 
DM particles are produced through t-channel annihilation process from constituent quarks 
of protons in the beam and the target nucleus. In case of the E613 we have a 400 GeV 
proton beam striking a tungsten target. The DM particles are produced from the process. 
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PP XX + The cross section for this hard process is calculated by incorporating the 
Lagrangian in eq.Q into MADGRAPH 5 [18] using FEYNRULES [T9|. The number of y’s 
produced is then given by 


dN _ r daipp^xx) 

dEde ~ * dEde 


( 2 ) 


where nt is the number density of nucleons inside the tungsten target, Nbeam is the number 
of protons in the beam that are incident on the tungsten target and Lt is the length of the 
tungsten target. The geometry of the E613 detector is such that only those DM particles 
are accepted for which the scattering angle 9 < 0.0134 |T5|. This is a conservative limit on 
the scattering angle compared to the detector acceptance in the original experiment where 
9 < 0.037. Thus we integrate eq.(|^ for 0 < 0 < 0.0134 to obtain DM distribution dN/dE. 


III. DEEP INELASTIC SCATTERING OF DIPOLAR DM 

Following [15] we place an energy cut of 20 GeV as the minimum energy required to 
register an event in the detector. Now in order to hnd the number of events we need to 
hnd scattering cross section of the DM produced in the experiment with the detector nuclei, 
having energy E^ where E^ > 20 GeV. Since these DM particles are produced in a collision 
with GM energy of 400 GeV, their scattering off the lead nuclei in the detector is deeply 
inelastic. The energetic DM particle therefore undergoes a deep inelastic scattering with 
the lead nucleus inside the detector via photon exchange. The DM couples to the photon 
through its electric/magnetic dipole moment. The initial DM momentum before scattering 
being k while that after scattering being k' the momentum transfer carried by the photon 
is q = k — k'. Now using the formalism of deep inelastic scattering of leptons we dehne the 

Bjorken scaling variable x = -, with = —q^ and v being the energy of the photon 

2m 

in the rest frame of the nucleus. With this the differential scattering cross section in terms 
of leptonic matrix element and hadronic matrix element is given by 

da ^ e^glpoie s 

dudQ'^ 167imN(E^ — m^) 

where Qdipoie = P, with /i being the magnetic dipole moment of DM and for DM interacting 
via the electric dipole moment, Qdipoie = E>. 

Now the leptonic current for the magnetic dipole moment interaction of DM corresponding 
to the Lagrangian in eq.([^ is 

Li^u = g 2 _ 2 _ 4^2 (g 2 ^/.v ^ ( 4 ) 

And similarity for the electric dipole moment interaction we have 

[4P;t- _ 2 + q^q^] (5) 

The hadronic matrix element Btyi, can be written in terms of structure functions such 
that we separate out contributions from longitudinally polarized photons and transversely 
polarized photons as [I^ 
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9fJ.u 2 

qz 


+ 2xa^^ ) Ft{x, + a^^FL{x, Q'^) 


where 




p ■ q + 2xmy 


Pm 


p-q 




Pm 


g 2 


-?M 


( 6 ) 


( 7 ) 


At the lowest order in perturbation theory the structure function Fl = 0 while Ft = 

Contracting the leptonic and hadronic tensors in eqs. Q and ([^ we have for the magnetic 
dipole interaction 


da = 


e^/i^ dudQ'^ V {2E — i>y 


z /2 + g 2 


IhTT 

Similarity for the electric dipole interaction we have 


- + 4m^ 


da = 


Ihvr E"^ — m‘1 Q'^ 


QH^E-uf ^2 . . 


1 /2+ g2 


g 2 - w 


'^xfq/A{x,Q^ 


^Xfq/A{X,Q^) 


( 8 ) 


(9) 


J Q 


For the nuclear parton distribution functions we use those provided by Hirai etah m- 
The above expression when integrated over u and gives the cross section for scattering 
of DM with nucleon inside the target. The limits of integration are as follows: 


Ecnt <v < E -m^ (10) 

g2 < g2 < 4(fc2 - Ev) - g2 (ii) 

1 ^9 2m2z/2 

where Qj =- = with k = E^ — m . 

k‘^ — Ev + J {kP' — Ei>Y — 

In addition to the above limits one also has the upper limit x < 1 which translates to 
g2 < 2mN^- From the cross section so obtained we can now write the mean free path of 
the DM particle as 


^ xN) 

where p is the number density of the nucleon inside the target material and a is the 
scattering cross section of DM from a nucleon. The probability of a DM particle scattering 
inside the detector is then given by P = 1 — and for a DM particle that behaves like a 
WIMP interacting through a weak dipole moment we can use the approximation P ~ L/A. 
Finally the number of expected events in the detector is 

Wv = J dE [Ppb(l - Ppe)] ^ (13) 

where Ppb is the probability of DM scattering inside the lead detector and Ppe is that inside 
the iron shielding. 
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FIG. 1: Left panel shows allowed parameter space (unshaded region) for electric dipole moment 
of DM in the — T) plane while right panel shows allowed parameter space for magnetic dipole 
moment of DM in the — 4 plane. 


IV. CONSTRAINTS ON DIPOLAR DM FROM BEAM DUMP (E613) AND 
OTHER EXPERIMENTS 


From the number of events calculated using eq.(13), we constrain the magnetic and 
electric dipole moments of DM in the DM mass range of 1 — 10 GeV. We use the interpretation 
in IZH also used in [15] of the experimental data in [161 IE] • We allow only those values of 
and 4 , V for which the number of expected events is < 180. The results for electric and 
magnetic dipole interactions of DM are shown in Fig. We have also added the bound from 
L3 experiment at LEP on the magnetic moment of DM [22]. The bounds on electric and 
magnetic dipole moments of DM from the analysis of L3 experiment are < 1.8 x 10“^^ e-cm 
or equivalently < 3.3 x 10“®4s. We see from the hgure that for electric dipole interaction of 
DM the allowed DM mass ranges between 1 — 3 GeV, however for the magnetic dipole case 
the DM mass lies between 1 — 8 GeV. In addition to the bound from L3 collaboration at LEP, 
recently constraints from solar physics on dipolar DM and similar models have also been 
studied [23l[26]. In [26] the bound on magnetic dipole moment of DM from helioseismological 
data is estimated to be 1.6 x 10“^^ e-cm, for DM mass < 4.3 GeV, which is quite similar to 
the bound from L3 collaboration. Also for momentum and/or velocity dependent scattering 
of DM studied in [23H25] which is relevant for dipolar DM model since it also has momentum 
and velocity dependence, the most favorable DM mass is found to be ~ 3 GeV. Thus we 
see that the constraints from beam dump experiments are in broad agreement with those 
from solar physics. Also the most stringent bounds from direct detection currently come 
from LUX collaboration [27], however this bound is considerably weak for low mass DM 
particularly in the 1 — 10 GeV range. As a result this bound is inconsequential for the 
results presented here. 


V. CONCLUSIONS 

We study the Dipolar DM model in the context of E613 beam dump experiment. We see 
that the constraint on electric and magnetic dipole interactions of DM from E613 experiment 
for light DM in the 1 — 10 GeV mass range is quite stringent. It is restricted in that mass range 
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mainly by the bound on magnetic dipole moment from the L3 collaboration at < 1.8 x 10“^^ 
e-cm or 3.3 x Bounds from solar physics data are also broadly compatible with 

this mass range and dipole moment. Thus the dipolar model of DM offers an alternative 
that is compatible with constraints from wide ranging experiments like beam dump or hxed 
target experiments as well as helioseismological data. In addition the low mass range of DM 
enables it to be compatible with the most stringent direct detection bounds from LUX [27]. 
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